The Journal of Experimental Biology

The ability of animals to adjust their heading, i.e. to turn, is essential for all walking animals. While several studies have addressed how leg movement or inter-leg coordination may change during turning, relatively little is known about how turning-related changes scale with turn magnitude. Here, we used spontaneous and visually induced turns of unrestrained walking stick insects to test (i) how high-level parameters of unrestrained turning scale with low-level parameters of leg movement, and (ii) the effect of visual guidance on turning parameters. To this end, we used a step change in stationary landmark position in an open-field arena to constrain timing and magnitude of target-directed turns. These visually guided turns were compared with spontaneous turns in an all-white condition. We show that visually induced turns were walked at a larger forward velocity and had fewer short steps than spontaneous turns. The scaling of turning responses was dominated by an increase in turning duration (factor 1.87) rather than turning speed (factor 1.32). Increased rotational velocity correlated with reduced forward velocity, though with flexible timing of both effects. These changes were accompanied by larger shifts in step direction, as well as an increased asymmetry of step types between inner and outer legs, suggesting a mix of distinct turning strategies, that depend on overall turn angle. Future models on six-legged locomotion should thus consider the incorporation of more than one mechanism to govern turning.

Predictions of skeletal muscle energy consumption under a diverse range of muscle contractile conditions are critical for improving our understanding of locomotion. Existing mathematical models, while capturing the mechanical dependence of energy consuming processes, neglect the time-dependent behaviour and recovery costs associated with regenerating ATP. This time-dependence is important for predicting the energetic response of muscles during repetitive or cyclical tasks like locomotion, where muscle undergoes many contraction cycles. This study presents a novel model to predict energetic rates based on physiological processes: Ca2+ transport costs, cross-bridge cycling costs, and ATP regeneration. Previous mathematical models include the dependence on Ca2+ transport and cross-bridge cycling, but neglect the time-dependent response and the subsequent recovery of ATP following the contraction. Model parameters were obtained from existing data on isolated muscle preparations, and predicted energetic rates were validated on separate datasets across a range of contractile conditions including dynamic, sub-maximal, and twitch contractions. The time-dependent model was able to capture the influence of contraction frequency on peak energetic rates and the time-course of energetic recovery observed experimentally. The model captures key physiological processes while maintaining a minimal number of free parameters and low computational cost. This enables generalisability across muscles and species, and implementation into larger scale musculoskeletal models.

Although muscle spindle sensory signals have been extensively studied, little is known about how and why muscle spindle firing is modulated by the central nervous system during movement. Specialized motor neurons to the muscle spindle, i.e. gamma motor neurons, can profoundly alter spindle firing during behavior, but technological limitations hinder our ability to record gamma motor and muscle spindle sensory signals during most behaviors. We used a biophysical model of a muscle spindle within a muscle-tendon unit to simulate how gamma drive may modulate muscle spindle Ia firing during locomotion. Based on a few available recordings from decerebrate animals, we demonstrate that our model, tuned to passive stretch conditions, can reproduce profound changes in muscle spindle firing in response to identical joint motions in locomotor vs. relaxed stretch conditions. Our model can discover phasic patterns of two types of gamma motor neuron drive based on recorded muscle spindle Ia firing and joint motion. By simulating perturbations, we conclude that: 1) sinusoidal activation of static gamma motor neurons during locomotion, encoding intended movement, modulates muscle spindle signals such that they act as sensorimotor feedback signals based on errors from the intended muscle fascicle length; 2) phasic on/off activation of dynamic gamma motor neurons during locomotion acts as an event detector, heightening muscle spindle Ia responses to discrete perturbations. As such, their muscle-within-muscle structure allows the muscle spindle to act as a highly tunable physical internal model of muscle state to guide movement. Our model supports proposed but as-yet-untested theories of muscle spindle function and offers a framework for extending the testing of muscle spindle function to active, behavioral conditions.

Parental care can have pervasive effects on offsprings neurodevelopment. Parent-offspring interactions are often modulated by the neuropeptide oxytocin, which is responsible for the development of social bonds. The development of the oxytocinergic system is dependent on the quality of parental care during the post-natal phase. However, it is yet unknown how post-natal direct interactions can influence the development of the oxytocinergic pathway. Here we tested how an obligate parental care behaviour, tadpole transport in poison frogs, influences the development of the oxytocinergic pathway. To this end, we quantified whole brain expression of oxytocin receptor and oxytocin precursor throughout three developmental stages of A. femoralis tadpoles, before, during and after tadpole transport. Our results show an overall downregulation during tadpole transport, which indicates that during transport tadpoles enter a dormant state to slow down development until they are placed in water. Interestingly, the expression of oxytocin precursor did not vary between the three developmental stages. This might indicate that oxytocin is being recruited during transport, but does not lead to neurodevelopmental changes. In sum, here we present the first evidence of a dormant state during tadpole transport which might be an adaptive response to the terrestrial reproduction in poison frogs.

The shift to terrestrial environments in ancestrally aquatic animals is often associated with key physiological and physical changes, including shifts in respiratory physiology and in some cases, even the evolution of completely novel respiratory structures. Examining how respiration operates across a gradient of submersion states in ancestrally aquatic terrestrial animals may shed light on how complex biological traits shift under different selective regimes. In this work, we begin exploring respiration in terrestrially-adapted land crabs that still use their gills to respire while underwater. We tested the relationship between aquatic respiratory rates, body size, and sex in red devil vampire crabs (Geosesarma hagen) at two ecologically-relevant temperatures. We found small females respire more than small males at 28{degrees}C, while large females respire more than large males at 21{degrees}C. Additionally, body size is a significant factor affecting respiratory rates of both sexes at 21{degrees}C and warmer temperatures significantly increase respiration in small crabs of both sexes. Interactions between these factors also led to emerging trends that can be explained by both physiological rules, such as reproductive investment and surface-to-volume ratios and heat transfer. We also report a temperature coefficient (Q10) of 1.52 for this species, showing an expected 52% change in respiratory and metabolic rate for every 10{degrees}C increase. This work also demonstrates the importance of understanding how and to what extent biological variables like sex and body size interact with abiotic environmental factors when measuring physiological traits in ectothermic invertebrate animals.

Animals fight for resources to obtain fitness benefits; most contests are intrasexual, and males tend to fight more than females. Although the genetic basis of male aggression is well studied, we know little about the genetic variation of female aggression. Female aggression varies with reproductive status and is potentially influenced not only by her genotype, but also by the genotype of her mate. Here we measured both male and female aggression in a set of Drosophila melanogaster inbred lines by competing each line against a standard competitor. Aggression varied among lines for both sexes, but male and female aggression were not correlated. Female aggression for many lines increased with mating, as expected, but not all lines changed aggression. However, when females were mated to males of different lines, male genotype did not affect the post-mating change in aggression, suggesting that ejaculate-mediated effects do not vary across these lines. The aggression level of the standard opponent was positively correlated with that of focal individuals indicating that individuals modulate their behavior according to the genotype of their opponent.

Muscle force sharing during locomotion is influenced by the mechanical demands of movement and the contractile properties of synergistic muscles. In cats, plantarflexor muscles exhibit distinct functional specialization, with the slow-fibred soleus maintaining relatively constant force across conditions while faster muscles such as the plantaris and gastrocnemius increase force production with increasing locomotor demand. However, it remains unclear whether similar force-sharing patterns occur in larger animals with different musculoskeletal designs. Therefore, the purpose of this study was to examine force sharing between the superficial digital flexor (SDF) and medial gastrocnemius (MG) muscles during treadmill locomotion in sheep. Tendon buckle force transducers were surgically implanted on the SDF and MG tendons of seven sheep, and in vivo muscle forces were recorded during walking and trotting across different speeds and inclines. Both muscles increased force with increasing speed and incline; however, speed had a substantially greater effect than incline. The SDF consistently produced greater absolute force than the MG across all conditions, whereas the MG exhibited slightly larger relative increases in force with increasing speed. Time to peak force decreased with increasing speed in both muscles, although the SDF reached peak force later in stance than the MG across conditions. In contrast to the distinct specialization observed in cats, neither muscle displayed a relatively condition-independent, soleus-like force contribution. These findings suggest that force sharing in sheep is more distributed across synergistic muscles and may reflect the influence of musculoskeletal design, tendon compliance, and mixed fibre-type composition on muscle function in larger species.

Rising ocean temperatures challenge ectothermic animals to maintain essential behaviors such as movement and feeding. We asked how a complete neuromuscular pathway preserves function when every component process responds differently to warming. In the pyloric system of the crab Cancer borealis, we simultaneously recorded motor nerve activity, muscle membrane potential, and contraction. Warming preserved rhythmic nerve activity and excitatory junctional potentials, but contraction declined and failed first. Fixed low-frequency stimulation, mimicking cold-temperature motor output, resulted in reduced contraction at warm temperatures, whereas higher-frequency stimulation, mimicking warm-temperature motor output, partially restored contraction. Warming hyperpolarized muscle fibers, moving them farther from contraction threshold, but also reduced input resistance, which together limited over-excitability. However, high-potassium stimulation revealed that the muscle contractile machinery remained functional. Thus, warming acts differently across levels, and overlapping compensatory mechanisms help preserve neuromuscular function across a wide range of temperatures. Significance statementCold-blooded animals that live in climates with significant seasonal changes in ambient temperature must have myriad mechanisms to function over a wide range of environmental conditions. We explore the effects of temperature at multiple levels of organization within the stomatogastric system of the crab, Cancer borealis. We find a series of compensatory mechanisms that cooperatively help maintain stable function despite the fact that the motor patterns, neuromuscular junctions and muscle functions are all differently temperature dependent.

Anoxia is a major stress for most vertebrates and frequently accompanies harsh winter conditions, particularly in species that spend much of the season frozen solid. North American freeze-tolerant wood frogs (Rana sylvatica) can survive several months without oxygen and endure whole-body freezing for up to eight months of the year, with [~]70% of total body water frozen as extracellular ice, yet revive when temperatures rise in spring. Survival depends on multiple adaptations, including tolerance of prolonged oxygen deprivation while frozen, when breathing and circulation are halted. A key strategy involves hepatic glycogen mobilization, producing large amounts of glucose that are distributed to tissues where it functions both as a cryoprotectant and as a substrate for anaerobic ATP production. The present study examines the role of histone lysine methylation and demethylation in regulating liver proteins under anoxic conditions. Relative protein expression of seven histone methyltransferases (ASH2L-S, ASH2L-L, RBBP5, SETD8, SMYD2, ESET, SETD1), six lysine demethylases (KDM1A, KDM3B, KDM4A, KDM4B, KDM5A, KDM5C), and eight histone marks (H3K4me1, H3K4me2, H3K9me3, H3K27me3, H3K36me3, H3K79me3, H4K20me1, H4K20me3) were evaluated in wood frog liver under control, 4-hour, and 24-hour anoxia exposures. The data indicate that histone lysine methylation and demethylation contribute significantly to transcriptional regulation under anoxia. Specifically, H3K4, H3K36, and H3K79 methylation were associated with transcriptional activation, whereas H3K9, H3K27, and H4K20 methylation correlated with transcriptional repression. These findings highlight the dynamic role of epigenetic regulation in supporting hypometabolism and stress adaptation in freeze-tolerant wood frogs.

Crown-of-thorns sea star (CoTS), Acanthaster cf. solaris, outbreaks are a major cause of hard coral cover decline across the west Pacific, threatening coral reefs. Coral taxa vary in susceptibility to CoTS predation from preferred (Acropora spp.) to non-preferred (Porites spp.), yet the mechanisms underlying these differences are poorly understood. We investigated coral defenses during an ongoing CoTS outbreak in Mo'orea, French Polynesia by examining gene expression (including putative toxin genes) in healthy and actively predated colonies of a preferred (Acropora hyacinthus) and a non-preferred (Porites sp.) coral prey species. During predation, A. hyacinthus exhibited molecular signatures of cellular stress responses involving oxidative stress signalling, inflammation, and tissue proteolysis. In contrast, Porites sp. showed enrichment of genes involved in mitochondrial metabolic adjustment and aerobic metabolism, suggesting metabolic compensation to maintain cellular function. Furthermore, A. hyacinthus demonstrated a reactive defense behaviour by differentially expressing toxins (e.g., kunitz-type neurotoxins) while Porites sp. employed constitutive expression of all putative toxins regardless of active predation, suggesting a proactive defense strategy. Together, these findings suggest that preferred and non-preferred coral prey exhibit fundamentally different molecular and defensive strategies during CoTS predation, shedding light on the evolutionary arms race between corals and their predators.

In a wide variety of animals, developmental crowding results in adults with smaller bodies. The crowding effect on body size in Drosophila melanogaster is canonically attributed to heightened competition for nutrition. However, whether other consequences of crowding also contribute to its effect on size remains an open question. We tested the relative contributions of nutritional competition, oxygen availability, and larval-generated metabolites to the crowding effect on size. We found that while nutrition explains most of the variation in body size due to crowding, oxygen also contributes in a sex- and nutrition-dependent manner. We found no evidence that larval-generated chemicals affect body size. These data confirm a widely suspected but untested role of nutrition in producing the crowding effect on size in D. melanogaster, while revealing an unexpected role of oxygen, and raise the possibility that behavior may be a mediator of density-dependent plasticity. Research HighlightsWe found that both nutrition and oxygen mediate the crowding effect on size in Drosophila melanogaster.

The age at which Atlantic salmon reaches sexual maturity shows a strong hereditary component associated with the vgll3a locus. The role of Vgll3 in maturation has remained unknown in vertebrates until recently, when it has been linked to pleiotropic roles in killifish, both delaying male maturation and affecting lifespan by protecting against cancer. As Atlantic salmon has two vgll3 paralogs, where only vgll3a has been associated with sexual maturation, it may provide a suitable model for studying the maturation-specific function of vgll3, as the other paralog may buffer for pleiotropic roles of vgll3. To address this, we used CRISPR/Cas9 to generate fish highly mutated in the vgll3a gene. We monitored their maturation and crossed highly mutated crispants to generate two year-classes of complete loss-of-function. All groups were reared under environmental conditions triggering early maturation in one-year-old males. We found a clear difference in the proportion of sexually maturing or mature fish between the different genotypes: in all experiments significantly fewer vgll3a-/- males entered puberty and reached final maturation compared to vgll3a+/- and vgll3a+/+ males. Furthermore, loss of vgll3a resulted in lower frequencies of maturation also in females. We conclude that Vgll3a stimulates maturation and that its complete removal significantly reduced maturation rates in both sexes in Atlantic salmon. Our findings also identify vgll3a as the causative gene in the locus associated with age at sexual maturity. Together, our findings support a new role for Vgll3 in initiating puberty in vertebrates and identifying salmon as a promising model for functional studies regarding the timing of sexual maturation.

Many insects enter diapause, a programmed state of developmental arrest that enables survival under adverse environmental conditions. In Drosophila melanogaster Meigen, 1830, diapause is characterized by reproductive arrest and reduced intestinal growth, accompanied by suppressed intestinal stem cell (ISC) activity. Juvenile Hormone (JH) promotes ISC proliferation under favorable conditions, but its capacity to modulate stem cell dynamics during cold-induced diapause remains unclear. Here, we investigated whether JH signaling can reactivate midgut remodeling in adult females maintained at 11. At this temperature, flies exhibited pronounced gut atrophy and elevated Phospho-histone H3 (PH3+) cell abundance, consistent with temperature-dependent G2/M phase arrest JH treatment significantly increased the proportion of Delta-positive progenitor cells in the anterior (R2) and posterior (R5) midgut regions at both 11 and 25, demonstrating that JH acts as a conserved mitogen for the ISC pool irrespective of thermal environment. A trend toward reduced PH3+ accumulation in the posterior midgut following JH treatment (p = 0.061) suggests possible facilitation of mitotic exit, though this effect did not reach statistical significance. Despite cellular-level changes, JH treatment did not restore overall gut size, indicating that the 72-84 hour exposure window was insufficient for subsequent tissue hypertrophy. Additionally, we identified a recurrent cold-induced pathology of gut distension, provisionally termed Lumen Obstruction Syndrome (LOS), which was independent of JH signaling. These findings reveal an uncoupling of JH-driven stem cell expansion from gross organ growth under diapause conditions, highlighting the selective sensitivity of the ISC compartment to endocrine signaling during environmental stress.

Despite extensive research on vestibular function in microgravity, particularly during orbital and parabolic flight exposure, several gaps remain regarding the spontaneous behavior of vestibular organs under non-terrestrial gravitoinertial conditions. In particular, semicircular canal function, typically assessed through vestibulo-ocular reflex (VOR) recordings, has yielded inconsistent findings, with reports describing either no effect or reduced performance in microgravity. Moreover, many of these studies are limited by methodological constraints that reduce the interpretability of their conclusions. To clarify these discrepancies, we evaluated horizontal and vertical VOR responses during parabolic flights to assess semicircular canal function under transient weightlessness. Participants were passively rotated at a constant frequency and amplitude during normogravity and microgravity phases, centered along the head vertical or inter-aural axis. Eye movements were recorded binocularly using infrared eye-tracking in darkness to eliminate visual influences, while participants were tightly restrained to minimize proprioceptive variability. Results show a reduction in VOR gain during microgravity in both axes, despite consistent rotational stimulation across gravity conditions. In addition, VOR gain remained reduced after parabolas in the horizontal plane, whereas vertical VOR performance was preserved. These are the first results to demonstrate an immediate alteration of semicircular canal function in weightlessness. Possible sources of the reduction in VOR performance in 0g are discussed. We also propose that the observed post-flight effects reflect a down-weighting of semicircular canal inputs during multisensory integration.

Many arthropods (insects and crustaceans) rely on their antennular chemosensory system to detect key environmental resources like food. While odor mediated food search is well studied in insects, characterization of crustacean chemosensory behavior has been limited by the long lifespans and large size of traditional crustacean model species. Here, we report the first characterizations of the food search behaviors of the genetically tractable amphipod crustacean, Parhyale hawaiensis. We find that Parhyale can locate an odorous food pellet, and predominantly approach food using direct, targeted swims from the arena walls. Removal of both first and second antennae dramatically reduced foraging success and impaired Parhyales ability to control take-off angle and maintain a stable heading during swims. Removal of the first or second antenna alone did not significantly disrupt foraging, and resulted in mild disruption of orientation phenotypes. Intact animals performed sharp turns near the location of the food pellet, which were observed when either first or second antenna were present, but not when all antennae were removed. Turns were longer and had higher average angular velocities following removal of either set of antennae, with full antenna removals representing the most extreme phenotype. In contrast with the long-held theory that the crustacean second antennae exclusively mediate contact chemosensation, we report that first- and second- antennae both contribute similarly to food localization and stabilization of locomotion in Parhyale in our behavioral paradigm. This work establishes Parhyale as an accessible model for studying olfactory behaviors in an aquatic arthropod.

The osteocyte network in bone is believed to play an important role for how bone tissues sense and respond to mechanical stimulation. Yet, bone adaptation to mechanical loads is often conceptualised as a simple response to mechanical stimuli, such as Wolffs law, which is based on mechanical variables only and takes no account of the cellular basis of mechanosensation. Wolffs law presumes the existence of a reference mechanical stimulus, the mechanical setpoint, above which bone is consolidated, and under which bone is removed. In this paper, we develop a theory of bone tissue sensing and adaptation based on osteocytes to provide new understanding of the role played by osteocyte signals in mechanical adaptation. In this theory, the mechanical setpoint of Frosts mechanostat is explicitly embodied as osteocyte properties involved in mechanotransduction. The mechanical setpoint is allowed to adapt due to the replacement of osteocytes during remodelling, making the setpoint space and time dependent. We propose a mathematical model to implement this new theory of bone adapation and present numerical simulations of this model to explore how mechanobiological response curves (effective Wolffs laws) are modulated by setpoint adaptation during remodelling. By accounting for varying osteocyte populations within bone tissue, we explore bone adaptation under osteocyte disruptions, which is particularly relevant to age-related bone loss. Our model suggests that biological disruptions of remodelling balance cannot always be compensated by mechanical feedback, and that setpoint adaptation during remodelling may have significant observable consequences, such as hysteresis in bone response signatures that resemble lazy zones.

In some reptiles, genetic and environmental sex determination interact whereby extreme incubation temperatures override genetic sex determination (GSD) to produce sex-reversed individuals. In one lizard with temperature-influenced GSD, the central bearded dragon, intron retention in the histone-modifier genes Kdm6b and Jarid2 has been implicated as a candidate signal linking temperature to sex. Equivalent intron retention is also present in two species with temperature-dependent sex determination, the red-eared slider turtle and the American alligator. The eastern three-lined skink, Bassiana duperreyi, represents another lizard with temperature induced sex reversal. It has an XY sex determination system in which low temperature incubation causes sex reversal of XX embryos to produce phenotypic males. In this study, we performed splice-aware analysis of RNA sequencing from hatchling brains of the three-lined skink. We investigated differences in alternative splicing and gene expression between the three sex conditions: XY males (XYm), XX females (XXf), and sex-reversed XX males (XXm). Sex reversal specific intron retention was observed in the gene, Ttll7, which only occurred in XXm and not in XYm or XXf. Intron retention in Ttll7 could alter the function of the encoded protein, a tubulin polyglutamylase, but its effect on sex reversal here is unknown. In addition, intron retention in the histone-modifier genes Jarid2 and Kdm6b occurred in all conditions. The presence of Kdm6b and Jarid2 intron retention in all sex conditions suggests that the pattern of intron retention in sex reversal in the eastern three-lined skink is distinct compared to the bearded dragon. We conclude that a different molecular pathway for sex reversal is induced in the three-lined skink, the details of which remain elusive.

Long-term winter dormancy in ectotherms (brumation) defines the annual cycle of many temperate-zone reptiles, yet the transcriptional regulation that supports survival across months of cold and aphagy remains poorly understood. We generated time-resolved transcriptomic profiles of liver and testis from male red-sided garter snakes (Thamnophis sirtalis parietalis) at five timepoints spanning the eight-month brumation cycle: pre-brumation, early, mid-, and late brumation, and post-arousal under continued aphagy. Time-course negative-binomial regression (maSigPro) followed by gene-set enrichment analysis identified 3,715 transcripts in liver and 5,828 in testis with significant temporal expression structure organized into five overarching temporal patterns: sustained downregulation, downregulation with post-arousal recovery, sustained upregulation, brumation-specific upregulation and cyclic modulation. Liver showed coordinated upregulation of fatty acid mobilization enzymes (ATGL, FOXO1, PPAR, CPT1A) and gluconeogenic regulators (CREBBP, PCK1) coincident with sustained low temperatures. Additionally, low temperature transcriptional activity was suggestive of a shift toward hepatic lipid mobilization and alanine-supported gluconeogenesis. Testis showed sustained suppression of meiosis, reproduction, and DNA-metabolism gene sets that did not fully recover at arousal consistent with this species dissociated reproductive pattern. Both tissues showed coordinated upregulation of stress-response pathways involving heat-shock proteins, HIF1 and a glutathione-based antioxidant defense. Interestingly, three vitellogenin transcripts and 17{beta}-hydroxysteroid dehydrogenases associated with estradiol-favoring steroid metabolism were upregulated in male liver during late brumation, which is not expected during natural physiology in adult males. Together these data support a framework in which temperature- and starvation-associated transcriptional programs contribute to survival of one of the longest, coldest brumations documented in a squamate. Summary statementA time-resolved transcriptomic analysis of liver and testis spanning eight months of winter brumation in Thamnophis sirtalis parietalis reveals gene expression patterns consistent with a temperature-associated shift toward hepatic lipid mobilization, sustained reproductive suppression, and vitellogenin response in males.

Computational Hill-type muscle models are widely used to simulate muscle force production because of their efficiency and physiological interpretability. However, their formulation relies on limiting assumptions, including debated multiscale simplifications, a simplified excitation-activation dynamics and an inability to capture slow and fast fibres. Moreover, existing Hill-type models remain insufficiently validated across physiological scales, fibre types, and contraction modes. We addressed these limitations by developing a multiscale fibre-type specific Hill-type neuromuscular actuator with mechanistic excitation-activation dynamics and systematically validated it against comprehensive experimental benchmarks. The model built upon a previously proposed motoneuron-driven actuator incorporating calcium-kinetics-based activation dynamics. The excitation-activation formulation was further refined to strengthen its physiological basis, while the contraction dynamics was extended by including an activation- and length-dependent force-velocity relationship, elastic tendon, passive elastic element, and the fibre-type-specific effects of yielding and sag. Validation was performed against four benchmark datasets spanning motor-unit and whole-muscle scales, including slow and fast fibres under both isometric and dynamic conditions. Experimental force traces were obtained from six muscles of rats and cats using a broad range of stimulation frequencies, muscle lengths, and imposed length changes, combining previous literature datasets with experiments performed ad hoc for this study. Overall, the model reproduced forces across all benchmark conditions, with mean absolute errors typically below 15% of the maximum isometric force, although larger errors were observed in specific submaximal and dynamic trials. The inclusion of physiologically based excitation-activation dynamics, together with yielding and sag, improved model performance under submaximal activation conditions. This study presents the first systematic validation of a single multiscale Hill-type neuromuscular actuator against comprehensive experimental motor unit and muscle force data, providing a benchmark framework for the development and assessment of future models. Author summarySkeletal muscles generate force through a complex sequence of events that links neural signals to muscle contraction. Because direct measurements are difficult to obtain, researchers often rely on computer models to investigate neuromuscular function and estimate muscle forces. However, most modelling approaches rely on simplifying assumptions about how force is generated across different biological scales, how muscles are activated, and how slow and fast muscle fibres behave. Moreover, they have not been validated against comprehensive experimental data. As a result, it remains unclear how accurately these models can reproduce muscle force across different physiological conditions. In this study, we established the first comprehensive set of experimental benchmarks spanning both motor-unit and whole-muscle scales, including slow and fast muscles under isometric and dynamic conditions. We used these benchmarks to validate a newly developed multiscale muscle model that explicitly represents the physiological pathway from neural stimulation to force production. The model incorporates experimentally based descriptions of calcium dynamics, activation, tendon elasticity, and fibre-type-specific contractile properties. We then compared simulated and experimental force responses across a wide range of stimulation frequencies, muscle lengths, and length-change conditions.

Ampullary electroreceptors are widespread across aquatic vertebrates. The purpose of sensing exogeneous electric fields is conserved across species but the implementations differ and the encoding mechanisms remain incompletely understood. We compared baseline and stimulus-driven response properties of ampullary electroreceptor afferents in the weakly electric fish Apteronotus leptorhynchus and Eigenmannia virescens. We find that their activity is well captured by an extended leaky integrate-and-fire model that generalizes across both species. The model shares similarities to a previous model of the tuberous electroreceptor afferents but further incorporates a low-pass pre-filtering and additional noise sources to reproduce the observed spectral response characteristics. The low-pass is essential to shape stimulus encoding in the high-frequency range. Accurate prediction of low-frequency stimulus encoding further requires two distinct noise sources: stimulus-independent white current noise and activity-dependent noise in the adaptation current, which is shaped by the adaptation time constant to yield effective pink noise dynamics. Using simulation-based inference, we trained a neural network to map model parameters to neuronal response features. This approach enables the generation of heterogeneous, biologically plausible model populations that may serve as a realistic input layer for studying neuronal processing on the next level. With this, we provide a unified and mechanistic model of ampullary electroreceptor encoding in these species and possibly beyond.